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Abstract 

The low-temperature thermal conductivity ko /T of d-wave superconductors is generally thought to attain a "universal" 
value independent of disorder at sufficiently low temperatures, providing an important measure of the magnitude of 
the gap slope near its nodes. We discuss situations in which this inference can break down because of competing 
order, and quasiparticle localization. Specifically, we study an inhomogeneous BCS mean field model with electronic 
correlations included via a Hartree approximation for the Hubbard interaction, and show that the suppression of 
Ko/T by localization effects can be strongly enhanced by magnetic moment formation around potential scatterers. 
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1. Introduction 

Thermal conductivity measurements at low tem- 
peratures T in the superconducting state have 
played an important role in strengthening the case 
for a d-w&vc BCS-likc description of quasiparticles 
in optimally doped cuprate superconductors. They 
are bulk probes of the superconducting state, unlike 
ARPES and STM, and can currently be performed 
at lower T than microwave experiments. One draw- 
back is the need to separate phonon and electron 
contributions, but in the cuprate superconductors 
an asymptotic linear term, kq/T ~ const which 
can be attributed solely to quasiparticles is almost 
always present at the lowest temperatures. After 
theoretical predictions of the universality of low-T 
quasiparticle transport in nodal superconductors [1] , 
experimental confirmation was obtained in a num- 
ber of optimally doped materials [2]. According to 
the theory, which relies on the disorder-averaged 
self-consistent T-matrix approximation (SCTMA), 
the low-T thermal conductivity is given by 
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where va is the gap slope at the node, and vf the 
Fermi velocity. This result has been shown to be 
insensitive to vertex corrections due to anisotropic 
scattering [3]. Since vp is generally well-known, 
Eq.(l) has been used to extract the gap slope for 
a number of cuprates at lower doping as well[4], 
leading to the conclusion that va increases with 
underdoping[4]. This conclusion is in apparent con- 
tradiction to recent ARPES experiments [5], so it 
is worthwhile to examine physical effects outside 
the framework of the SCTMA which could lead to 
a suppression of kq/T and thereby to the possible 
erroneous conclusion based on Eq. (1) that va was 
increasing. 

There are several effects known to lead to a sup- 
pression of Ko/T. Localization effects were discussed 
in this context in Ref. [6] , and effects of bulk subdom- 
inant competing orders have been shown to suppress 
kq/T but do not immediately eliminate it, despite 
the removal of the d-wave nodes[7]. Here, we inves- 
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tigate the effects on k(T) by local impurity-induced 
moments relevant e.g. to the underdoped regime of 
La 2 -xSr x Cu04 (LSCO)[9]. 

2. Magnetic correlations. 

The model used to study disordered d-wave su- 
perconductors with magnetic correlations is: 

H = - 1 ]T clc ja + £ (& l3 c\Al + H ' c -) 

(ij)<r (ij) 

- ^2 ^hia + U ^2 flia{flia) + Vj,n irT . (2) 
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Here, c ia creates an electron of spin a on the site i 
on a two-dimensional lattice, /i is the chemical po- 
tential and fiicr — c irr Cia is the particle number op- 
erator. In general, magnetic order induced by U will 
compete with the superconducting order Ay and 
lead to a bulk magnetic state above some critical 
Uco- However, even for U c < U < U c q, disorder V{ ^ 
caused by point-like impurities can induce local- 
ized S = 1/2 states with staggered magnetization. 
Here, U c is the critical U necessary for local moment 
formation, and will in general depend on the impu- 
rity strength. In this paper we study impurity con- 
figurations with 2% randomly positioned point-like 
impurities of strength Vj, = 35< which is in the uni- 
tary limit. We have solved Eq. (2) selfconsistently 
on 40 x 40 lattices with hole doping 6 = 12%, and 
calculated the electronic k(T) as outlined in Ref. [6] . 

In Fig. 1, we show the thermal conductivity aver- 
aged over 20 different random impurity configura- 
tions. For [7 = the results agree well with those 
obtained previously for the dilute impurity limit [6]: 
k(T)/T = Koo + aT 2 , and the spatial inhomogeneity 
of Ay only slightly modifies a as seen by comparing 
the SC and NSC results. For the bulk system we have 
Ay = OAt on each link resulting in the universal 
value Koo = 1/3 (vf/va + va/vf) — 0.967, which is 
nicely reproduced by the numerics. For < U < U c , 
\/i : (riii) — (riii) ano - ^ enters as a chemical poten- 
tial shift which is compensated by a modified /i (to 
get S = 12%) and consequently k(T) is identical to 
the C/ = result. In the regime U c < U < C/ c0 , local 
moments are formed and the corresponding thermal 
conductivity is also shown in Fig. 1. As seen, when U 
increases, n(T)/T (and Kq/T) is continuously sup- 
pressed, and kq/T will eventually vanish in the bulk 
insulating state around U ~ U c q. The origin of the 
suppressed thermal conductivity can be traced to a 
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Fig. 1. k(T) /Tkqo versus T for 2% point-like impurities with 
U = (non-selfconsistent (NSC) (□), selfconsistent (SC) 
(*)), U = 3.2t (A), and U = 3At (o). The finite system size 
restricts us to study the region T/t > 0.02. Inset: Selfcon- 
sistent DOS for the cases U = 0.0 (*) and V = 3.2t (A). 

reduction of the low-energy density of states (DOS) 
by the local magnetic moments as shown in the inset 
in Fig. 1. Thus, impurity-driven local moment for- 
mation in the underdoped regime may explain the 
doping dependence of kq/T measured in LSCO[10]. 
In addition, since the low-T thermal conductivity is 
severely suppressed, measured values of koo an d a 
naive use of the "universal" clean <i-wave result in 
Eq. (1) would lead to an erroneous estimate the su- 
perconducting gap. 
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